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Abstract 

Using a high statistics sample of photo-produced charm particles from the FOCUS 
experiment at Fermilab, we report results of a search for eight rare and Standard- 
Model-forbidden decays: Df — > (with h = Tr,K). Improvement over 

previous results by a factor of 1.7-14 is realized. Our branching ratio upper limit 
D+ 7r+/i-^+ of 8.8 X 10-6 at the 90% C.L. is below the current MSSM R-Parity 
violating constraint. 



1 Introduction 

The search for rare and forbidden decays of charm particles is enticing since 
Standard Model (SM) predictions for interesting decays tend to be beyond 
the reach of current experiments, and a signal is an indication of unexpected 
physics. 

Standard Model predictions [1,2,3] for the branching ratios of rare decays 
D^,Df — >• hf^ ii~ 11'^ (with h = n,K) are dominated by long range effects 
which are notoriously difficult to calculate. Even so, the differences between 
the three predictions in [1,2,3], although using individual treatments for the 

* See http://www-focus.fnal.gov/authors.html for additional author information. 
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decay spectra, are quite manageable. For example, the predicted integrated 
rate for n^fi^fi^ varies by only a factor of 2 while the experimental 

limits are a factor of 5-10 away. Minimal Supersymmetric Standard Model 
(MSSM) R-Parity violating extensions can significantly increase this rate. Ex- 
perimental results for D'^ — > tt"^ currently set the best constraint for 
the product of the MSSM R-Parity violating couplings: |A^2fe^2ifer < 0.004 [1]. 
Until experiments reach the SM limit for these rare decays, a signal indicates 
new physics or a needed refinement in the interpretation of the SM. Decays of 
the form D+, — > (with h — 7r,K) are forbidden in the SM since 

they violate lepton number conservation, and a signal in these modes is direct 
evidence of new physics. 

In this paper, we present new upper limits for the branching ratios of 3-body 
di-muonic decays of the D"*" and mesons mentioned above. Unless specif- 
ically mentioned, all results include a lower limit of at the 90% C.L.. Our 
results represent a factor of 1.7-14 improvement over previous experimental 
limits [4,5]. The result for the branching ratio upper limits of n^fi^fi^ 
of 8.8 X 10-6 @ 90% C.L. and Df ir+n+fi- of 2.6 x lO^^ @ 90% C.L. are both 
within a factor of 5 of the Standard Model (long range) predictions 1.9 x 10"^ 
[1] and 6.1 X 10"^ [2] respectively. The result for the branching ratio upper 
hmit D+ of 8.8 x 10"^ @ 90% C.L. is below the MSSM R-Parity 

violating constraint [1]. 

The data for this analysis were collected using the Wideband photoproduction 
experiment FOCUS during the 1996-1997 fixed-target run at Fermilab. The 
FOCUS detector is a large aperture, fixed-target spectrometer with excellent 
vertexing and particle identification used to measure the interactions of high 
energy photons on a segmented BeO target. The FOCUS beamline [6] and 
detector [7,8,9,10] have been described elsewhere. 



2 Event Selection 

We look for D^s through the 3-body decay chains D'^, — > h^n^n'^ (where 
the h represents a pion or a kaon) for rare decays and K~ti^ti^ 
or — * K~ K^-K^ for normalization (charge conjugate modes are implied 
throughout this paper). In order to search for the set of cuts that provides 
signal optimization, we place initial (loose) requirements on the reconstructed 
data to produce a base sample, and then we place a series of (tighter) cuts on 
the base sample. The loose requirements consist of acceptance, momentum, 
vertexing, and particle identification cuts. Note that for all cuts, care is taken 
to ensure that the normalization modes receive the same cuts as the di-muon 
modes where possible. 
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Dtic to the finite lifetime and Lorentz boost of charm candidates, the primary 
interaction vertex and secondary decay vertex can have a significant separation 
along the beam direction. Secondary vertices are formed from 3 candidate 
tracks, and the resulting momentum vector of the 3 tracks is used as a seed 
to search for a primary vertex [7] . We require that the confidence level of the 
secondary vertex fit (DCL) exceed 1%, the confidence level of the primary 
vertex fit exceed 1%, the significance of separation between the primary and 
secondary vertices (^/cr^) exceed 5, and the confidence level that any of the 
secondary tracks is consistent with the primary vertex (ISOl) be less than 
10%. The latter cut is included to remove the contamination from D*+ decays 
and other tracks originating from the primary that could be confused with 
secondary tracks. 

We use the Cerenkov system [8] to identify pions and kaons. For each track, 
Wohs — —2 log (£) is computed, where C is the hkelihood that a track is 
consistent with a given particle hypothesis. For a candidate kaon, we require 
W^obs(7r) — Wohs{K) (kaonicity) be greater than 0.5. For a candidate pion in a 
rare mode, we require PVobs(-^) ~ W^obs(7'") (pionicity) be greater than -15. 

Muon candidates are required to be within the acceptance of either of the 2 
muon systems in FOCUS [10]. We require the tracks in a normalizing mode 
corresponding to muon tracks in a rare decay mode be in the acceptance of 
one of the muon systems as well. Since the rate of muon misidentification 
increases at low momentum, we place a requirement that the momentum of 
"muon" tracks within the acceptance of the wide angle (outer) muon system 
t'c greater than 6 GeV/c and those within the low angle (inner) muon 
system (p^""icr^ greater than 8 GeV/c. All muon candidates are required 
to have associated hits in either muon system sufficient to meet a minimum 
confidence level, /icl, of 1% for the muon hypothesis, and must pass additional 
muon cuts depending on the system traversed. For the outer muon system, 
muon candidates must traverse a minimum of 150 cm of material, produce hits 
in 2 of 3 planes of the detector, and these 2 (or more) hits (called a cluster) 
must not be shared by the other muon candidate. For the inner system, at least 
4 of 6 planes of the detector must record hits consistent with the candidate 
track, no more than 2 of these hits can be shared with the other candidate 
muon track, and a fit to the hypothesis that the inner muon candidate track 
had the same measured momentum in both magnets traversed was required to 
exceed 1%. This last cut is used to reduce contamination from particles that 
decay and produce a muon as they traverse the spectrometer. This cut is also 
applied to the pions in the normalization modes, and the lowest momentum 
kaon, when possible, for the normalization mode. Finally, all other tracks 
in the event are fit to the muon hypothesis using the hits from a candidate 
muon, and the highest confidence level from the fits, ISO^, is saved. No ISO^ 
cut was required for the base sample, but a 10% ISO^^ is included in the set 
of cuts used for sensitivity optimization. 
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For the normalization signal, D'^ K~^K~7r^, a cut was applied to reduce 
the reflection when one of the pions from K^ii^n^ is misidentified as 

a kaon. Under the hypothesis that the same sign kaon track is assumed to 
be pion, the invariant mass is calculated. If the new invariant mass is within 3 
standard deviations of the D'^ mass, the kaonicity of the same sign kaon must 
exceed 6. Additionally, we require that the reconstructed K'^K" invariant 
mass be within 3 standard deviations of the (j) mass. 

Our analysis methodology (Section 3) requires a base sample of events of suf- 
ficiently small size (150 events) for a reasonable processing time. Base samples 
were obtained by applying the minimum l/a^ cut in the range of 5 to 8 which 
brought the sample size below 150. The series of cuts applied to these sam- 
ples was arranged into a grid. This cut grid was based on kinematic variables, 
particle ID algorithm results, vertex quality, and event topology. For example, 
once the lower cut in f/cr^ was determined, the l/a^ cut was allowed to vary in 
steps of 2 units up to a maximum of 21 or 22. The cuts used have been shown 
to be effective for other charm decays besides those presented in this analysis. 

Since the kaonicity cut can be applied to the kaon modes identically, there is 
an inherent increase in the size of the cut grid when the pionicity cut is used 
for the pion modes. To keep the size of the cut grid roughly the same, and to 
prevent memory overflows in software that was difficult to modify, we reduced 
the size of the cut grid for the pion modes. The grid (see Table 1) includes 
cuts on £/(J^, ISOl, DCL, kaonicity, pionicity, ixcli ISO^, and the momentum 
of muon candidates. The normalizing modes used to compute the branching 
ratios for and are shown in Figure 1 for the loosest cuts in the grid. 




Fig. 1. Modes used to normalize the rare decay modes. The soUd histograms repre- 
sent the loosest cuts employed in the analysis, while the superimposed cross-hatched 
histograms represent data which has had the tightest cuts used in the analysis ap- 
plied. Notice the large reduction in background relative to the signal for both the 

and the modes over the range of cuts used. 

In order to perform signal optimizations, calculate efficiencies, and estimate 
yields, there is a distinction made between the signal region, where events for 
the desired mode are expected to occur, and the background sidebands, where 
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Table 1 

Analysis cuts used in the cut grid. Variables are described in the text. The best cut 
on average represents a point on the cut grid used in a systematic check of our result 
that is described in the "Systematic Checks and Results" section of the paper. Cuts 
indicated by { } are appHed only to the kaon modes to keep the cut grid about 
the same size for kaon and pion modes. Notice that the cuts removed for the pion 
modes are chosen either very far from the "Best Cut" (explained later in the text), 
or represent a small reduction in the stepping of a cut that varies logarithmically 
(ISOl). 



Variable 


Cut Values Used in the Grid 


Best Cut on Average 




> 5 ^ 22 


> 13 (Z?+), 10 (L>+) 


ISOl 


< 0.1, {0.03}, 0.01, {0.003}, 0.001 


< 0.1 


DCL 


> 1%,2%,4% 


> 1% (I?+), 2% (D+) 


Kaonicity (kaon modes) 


> 0.5,1.0,2.0 


> 1.0 


Pionicity (pion modes) 


> -15,-3,-1 


> -3 


^i'CL 


> 1%, 5%, 10% 


>5% 


iso^ 


< 0.10,1.0 


< 1.0 


P^'"^"^-^ (GeV/c) 


> 8, 9, 10, 11, 12, {14} 


> 9 


P^""*^"^ (GeV/c) 


> 6, 7, 8, 9, {10} 


> 7 



the amount of signal is expected to be minimal. The background sidebands 
are used to estimate the amount of contamination in the signal region. We 
define the signal region to be within 2 standard deviations of the nominal 
reconstructed parent particle mass (i.e. either the (± 20 McV/c^) or the 
(± 18 McV/c^)) and the background region to be any invariant mass 
reconstructed outside the signal regions between 1.7 and 2.1 GeV/c^. An ex- 
ception is made for D"*", /.i^/j'^. For K~fi~^/j^ the background 
sidebands are required to be between 1.75 and 2.1 GeV/c^. This approximately 
splits the expected contribution of i^~7r^7r^, where the 2 pions arc both 
mis-identified as muons, equally between signal and background regions. For 
Df — > K~ ^'^jx^ we require the lower sideband to begin 2 standard deviations 
above the mass. This effectively removes the — > K~t:^t:^ signal that 
comes from mis-identified muons. The centroid and standard deviation for 
each mode are determined by fitting the reconstructed parent particle mass 
in the normalization mode. Since there is a small shift in the mass centroid 
between data and Monte Carlo, we use the data to determine the regions for 
data and Monte Carlo to determine the regions for Monte Carlo. 

The shape of the background in each rare mode is determined using a large 
sample of photoproduction Monte Carlo events where all charm species and 
known decay modes are simulated. The shape of the background is used to 



6 



determine r, which is the ratio of the number of background events expected 
in the sideband region to the number of background events expected in the 
signal region. An average r is computed for each mode from all the r's in 
the cut grid for that particular mode. We find that the ratio of Monte Carlo 
efficiencies and r are stable over the cut combinations in the grid, but we 
require at least 10 surviving Monte Carlo events to determine each r used 
in the average to prevent variations due to low Monte Carlo statistics in the 
signal region. 

These definitions are used in the determination of the branching ratio limits 
in the analysis described below. 



3 Analysis 

The analysis technique emphasized a careful approach to the treatment of 
backgrounds in a limited statistics analysis. The "blind" approach, to select 
cuts that optimize signal efficiency relative to background sidebands, may 
still lead to a downward fluctuation of the sidebands relative to the masked off 
"signal" region and a more conservative limit on average [11]. Further, authors 
frequently use the technique outlined by Feldman and Cousins [12] to calculate 
the confidence levels used in the calculation of their limits. The Feldman- 
Cousins approach does not explicitly include fluctuations of the background. 
Indeed, the Particle Data Group [13] suggests presenting a measure of the 
experimental sensitivity, deflned as the average 90% confidence upper limit 
of an ensemble of experiments with the expected background and no true 
signal, in addition to the reported limit whenever experiments quote a result. 
None of the methods suggested in the PDG, including the Cousins-Highland 
method for including systematic errors in upper limits [14], properly deal with 
fluctuations in the background and bias in selecting the data. 

For this analysis, we chose a method suggested by Rolke and Lopez [15] which 
includes the background fluctuations directly into the calculation of the likeh- 
hood. The composite Poisson probability of flnding x events in a signal region 
and y events in background sidebands given a signal // and a background h is: 

P^,b{x,y)^ y\ 

where r is the expected ratio of the number of background events in the side- 
band regions to the signal region. Rolke and Lopez have shown that including 
the second Poisson factor in this expression leads to confidence intervals with 
better coverage than those of Feldman-Cousins who only consider the first 
factor. 
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In a second paper Rolke and Lopez have shown that bias can also be intro- 
duced during the selection of optimal cuts [11]. If a single cut is chosen based 
partly on the level of background in the sidebands (as is typical), there is a 
tendency to optimize on downward fluctuations and, hence, to underestimate 
the background level in the signal region. The resultant limits from such an 
"optimized" analysis, even though carried out in a "blind" fashion, will not 
have the correct coverage. 

In order to reduce the bias due to selection, Rolke and Lopez suggest the data 
be sampled using a "Dual Bootstrap" method. In a bootstrap, the experimen- 
tally observed data set of N events is used to create an ensemble of many 
different N-event experiments or data sets, obtained by random sampling of 
the original data set allowing repeated events. In the Dual Bootstrap, two 
independently bootstrapped data sets are created. One set is used to optimize 
the cuts which are then applied to the second set in order to calculate the 
confidence intervals. This procedure is repeated 10,000 times and the median 
value for the limits is the final result. The two bootstrap data sets are suf- 
ficiently independent that the background estimate from the second is very 
nearly unbiased. 

We use the experimental sensitivity (see Eqn. 2 and Eqn. 3) as our figure 
of merit to optimize cuts. A matrix or grid of possible data quality selec- 
tion cuts are applied to the first bootstrap data set. The point in the multi- 
dimensional cut grid which has the best (smallest) sensitivity is applied to the 
second bootstrap data set. Limits for the branching ratio and a new sensitivity 
are computed from this second set. The Dual Dootstrap procedure is shown 
schematically in Figure 2. 

For a given r with y sideband events, the average 90% confidence upper limit 
for the number of events in the signal region when there is no true signal, 
S'^(y), can be calculated from the Rolke- Lopez [15] 90% upper limit table as: 

oo 

Sr{y) = Y.Ur{x,y)-Py/r{x) (2) 

where Ur{x,y) is the upper limit of the signal, and Py/T{x) is the Poisson 
probability of x when the expected background is y/r. 

The experimental sensitivity is: 

S ( ^ 

Sensitivity = BRnormTT — e (3) 

-'norm 

where the branching ratio of the normalization mode, BRnorm, comes from the 
PDG [16], Fnorm IS the yield of the normalization mode (determined from a 
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sideband subtraction), and e is the ratio, determined using Monte Carlo, of the 
normahzation mode efficiency divided by the rare decay mode efficiency. For 
different bootstrap data sets and different cut sets, y, r, ynorm, Sriy), and e 
can be different. The sensitivity meets the requirement of a blind analysis, i.e. 
it does not depend on the number of events observed in the signal region. In a 
Dual Bootstrap procedure one sensitivity is calculated as the best sensitivity 
for the first data set and a second (less-biased) sensitivity is calculated when 
the "best" cuts are applied to the second data set. 

The 90% confidence upper limit for the rare branching ratio is: 

^rarc 

Upper Limit = BRnorm TT^e (4) 



where l^are = Ur{z,y) is the Rolke-Lopez 90% confidence upper limit for the 
signal yield given z events in the signal region. The lower limit has a similar 
expression. 



Rare Event Candidates 

100 events 



. Randomly sample 



Randomly sample 



1st Bootstrap Sample 



100 events 



2nd Bootstrap Sample 



100 events 




Save Cut Combination 



Use 

1st Bootstrap Cut 
on 

2nd Bootstrap Data 



Save Result 



Perforin Dual Bootstrap 
10,000 times 



10,000 Sensitivities 
& Branching Ratios 



Fig. 2. A flowchart for the Dual Bootstrap on a sample size of 100 events. Note the 
parallel structure of the sampling with replacement that separates the cut selection 
that optimizes sensitivity from the calculation of the sensitivities and branching 
ratios used to determine the final result. Thus, only the data sidebands and the 
expectation of the shape of the background is used in the cut selection, even though 
events may be chosen from throughout the data set during the bootstrap process. 
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The salient features of parameters coming from the Dual Bootstrap analysis 
are illustrated in Figure 3 for the decay mode — > 7i~^fi~fi~^. Several fea- 
tures are worth mentioning. The spread of the branching ratio upper limits is 
in general larger than the spread of the best sensitivities. This happens since 
the spread of the branching ratio upper limits is dominated by fluctuations 
in the signal region while the spread in the best sensitivities stems from the 
larger background region used in the estimate. This is a distinct advantage 
when determining r with large sidebands rather than smaller sidebands where 
the fluctuations in the background can become more problematic. The single 
bootstrap sensitivities sample the very lowest end of the sensitivity spectra, 
and fluctuations often place the best bootstrap sensitivity below any calculated 
from the data directly. Notice though, that the second bootstrap sensitivities 
remove the bias and give a median sensitivity somewhat above the minimum 
expected by the data. This is a safeguard against choosing a single cut that 
produces an outlier or poor estimation of the true sensitivity. Thus, the Dual 
Bootstrap method allows us to optimize the sensitivity for each decay mode 
while retaining correct coverage. Also notice that the Dual Bootstrap branch- 
ing ratio upper limits can cover a larger spread than the original branching 
ratio upper limits shown for all cut combinations. This can be understood if 
one realizes that for a small number of events in the signal region, the back- 
ground can vary considerably in the second bootstrap, and the confidence 
interval calculation amplifies the effect, producing a larger spread. 



xlO 
1400 - 



- Dual Bootstrap Median Sensitivity 



- Dual Bootstrap Median B.R. Upper Limit 



—All Cut Combinations Sensitivity 



- Single Bootstrap Best Sensitivities 



- Dual Bootstrap Sensitivities 




- Dual Bootstrap Median Sensitivity 



- Dual Bootstrap Median B.R. Upper Limit 



All Cut Combinations B.R. Upper Limits 



Dual Bootstrap B.R. Upper Limits 



"0 0.05 1 0.15 0.2 0.25 3 0.35 0A_^ "O 0.05 1 0.15 0.2 0.25 3 0.35 0.4_^ 

Sensitivity ^ 10 Branching Ratio ^ 10 

Fig. 3. The distribution of results used in the determination of the branching ratio 
confidence interval for the decay mode — > Tr^^~^~^. In plot a), the sensitivity 
distributions are shown, and in plot b) the branching ratio distributions are shown. 



10 



4 Systematic Checks and Results 



The largest sources of systematic uncertainty in this analysis are estimated 
to come from uncertainties in the Monte Carlo simulation and uncertainty in 
the branching ratios used for the normalizing modes. One source of systematic 
uncertainty is included specifically in the upper limits through the r param- 
eter, while the other sources of systematic uncertainty were included using 
the technique outhned in [14]. Using this method, the increase, AC/„, in the 
Poisson upper limit on the estimate for the number of rare decay events is: 

1 o 9 U^j +b — n 



where U^^ represents the Rolke-Lopez 90% confidence upper hmit for the 
mean number of signal events, l^are, b is the predicted background in the 
signal region, n is the number of events found in the signal region and the 
total relative systematic uncertainty is cTr- For example, uncertainty in the 
normalizing branching ratios and efficiencies used in the calculation of the rare 
branching ratios is translated to a percent or relative error in the estimation 
of Y 

^ rare ■ 

The relative systematic uncertainty from the normalizing branching ratios, 
cpDG) comes from the PDG [16]. The relative systematic uncertainty stemming 
from the simulation of the data comes from the simulation of the experimental 
trigger, (iTrigger, and the estimation of the efficiency of the outer muon system, 
(j^iD. Since the FOCUS trigger requires a minimum energy be deposited in 
the calorimctry, muons and hadrons will deposit very different energies, and 
the trigger simulation must account for any difference. The difference between 
three very different simulations is used to estimate this uncertainty: a full 
GEANT [17] simulation of the calorimetry, a pre-stored shower library gen- 
erated with GEANT which selects the calorimetry response based on particle 
types, energies, and locations, and a simulation based on the parameterized 
average response in data of the calorimetry based on incident particle types 
and energies. The relative systematic error due to the outer muon identifica- 
tion is estimated by looking at the difference in the Monte Carlo rare decay 
efficiency for 2 different estimations of the the outer muon identification effi- 
ciency One method employs the overlap between the inner and outer systems 
(very parallel muons coming from far upstream of the experiment can impact 
both systems) , and the other method uses 2 hits in the outer system to predict 
a third hit. The final source of systematic error is due to the uncertainty in 
the modelling of the muon misidentification used when r is calculated for the 
K^fi^fi^ decay mode. This uncertainty is estimated by boosting the 
contribution of — > K~7r~^7r~^ in the photoproduction Monte Carlo by twice 
the amount needed to match the amount of — > K'tt'^tt'^ seen when one 
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of the pions is misidentified by a muon. The more conservative r is then used. 
The small statistical errors from the ratio of Monte Carlo efficiencies and the 
error in the yield of the normalization modes did not contribute significantly 
to the systematic error. 

The sources of relative systematic error for each mode are shown in Table 2. 
The total relative systematic uncertainty is obtained by adding all the contri- 
butions in quadrature. The effect on the rare branching ratio is calculated for 
each bootstrap sample and is naturally included in the ensemble result. 

Table 2 

Contributions to the relative systematic uncertainty, o"r, in %. 



Decay Mode 


•^Trigger 




CPDG 




D+ K+fi+fi- 


2.8 


1.9 


6.7 


7.5 




2.7 


2.6 


6.7 


7.7 




2.5 


2.7 


6.7 


7.6 




2.0 


2.6 


6.7 


7.5 


Dt ^ K+^^+^^- 


3.0 


1.9 


27.3 


27.5 


D+ ^ K-^^+^^+ 


2.3 


2.5 


27.3 


27.5 


D+ ^ 7r+/x+//- 


3.6 


2.7 


27.3 


27.7 


D+ ^ 7r-//+/x+ 


3.0 


2.8 


27.3 


27.6 



To compare the Dual Bootstrap results to a more traditional "blind" analysis, 
another technique was used that selected a unique set of cuts, or point on 
the cut grid. Since the and Df lifetimes and production topology differ, 
a separate point on the cut grid was determined for each. The cuts used to 
determine the best sensitivities in the first bootstrap are saved for all four 
rare modes of a parent particle. A point in the multi-dimensional cut grid is 
determined by choosing cuts closest to the average value of each saved cut (see 
Figure 4). The cuts represented by these 2 cut grid points, one for the and 
one for the D^, are then applied to the respective modes once in the spirit of 
a more traditional "blind" analysis. A branching ratio limit is computed using 
the resultant data histogram and the previous definitions for the signal region, 
the background sidebands, and r. The best average cuts are shown in Table 1, 
and the data histograms resulting from this check are shown in Figures 5 and 
6. A comparison was also made between the confidence limit calculated using 
the Rolke- Lopez method [15] and the Feldman-Cousins method [12]. Little 
difference was seen. We stress that these checks are provided as a convenience 
to the reader. As stated previously, the methods of Rolke and Lopez [11,15] 
have been demonstrated to provide correct coverage, whereas the coverage of 
the checks mentioned has either not been studied or, in the case of Feldman- 
Cousins where background fluctuations are not considered, has been shown to 
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have incorrect coverage. 



Rare Event Candidates 
(Same Topology/Parent) 

100 events 



Randomly sample 



Bootstrap Sample 




100 events 


►< 

Apply Cut Grid 




Save Cut Combination 



Perform Bootstrap 
10,000 times 



Find Average Cut from Ensemble 
Choose Closest Point on Grid 



Apply Chosen Point on Cut Grid to Data 
Calculate Branching Ratio Upper Limit 



Fig. 4. A flowchart for the cut bootstrap on a sample size of 100 events. This hgure 
should be compared to Figure 2. In this technique we determine an average "best 
cut" for the data using only the data sidebands and the expected shape of the 
background. 



No significant evidence for the observation of any of the rare decay modes 
was seen. All modes except K~fi^fj,~^ had a 90% lower limit for the 

branching ratio of zero. For — >■ K^n^fi^, the 96% lower limit was zero. The 
results of the analysis are presented in Table 3 below. There is good agreement 
between the Dual Bootstrap branching ratios (for both the Feldman- Cousins 
and Rolke-Lopez limits), the sensitivities, and the single cut systematic checks. 



5 Summary and Conclusions 



The FOCUS results from this analysis are shown in comparison to previous 
best results and recent theory in Table 4. Our results are a substantial im- 
provement over previous results [4,5] and the FOCUS result for the branching 
ratio upper limit — »■ 7r+/x+/x~ of 8.8 x 10"^ @ 90% C.L. is lower than the 
current MSSM R-Parity violating constraint [1] for this mode. 
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Fig. 5. Data used in the single cut systematic check for the decay modes. Note 
that the r's shown on the plots are the same as those used for the Dual Bootstrap 
analysis. The solid histogram entries correspond to events in the signal region. 
The cross-hatched areas to either side of the normalization mode signal correspond 
to the data used for the sideband subtraction. The singly hatched areas in the 
di-muon mode histograms correspond to the signal region, while the cross-hatched 
areas correspond to an excluded region. All the other data and area shown in the 
di-muon histograms are used for the background estimate. 
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Fig. 6. Data used in the single cut systematic check for the decay modes. Note 
that the r's shown on the plots are the same as those used for the Dual Bootstrap 
analysis. The solid histogram entries correspond to events in the signal region. 
The cross-hatched areas to either side of the normalization mode signal correspond 
to the data used for the sideband subtraction. The singly hatched areas in the 
di-muon mode histograms correspond to the signal region, while the cross-hatched 
areas correspond to an excluded region. All the other data and area shown in the 
di-muon histograms are used for the background estimate. 
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Table 3 

FOCUS results with and without incorporated systematic errors for the modes 
shown. Each number represents a 90% confidence upper hmit for the brancing ra- 
tio of the decay mode listed. F-C represents the Feldman-Cousins 90% confidence 
upper hmit. R-L represents the Rolke-Lopez 90% confidence upper hmit. Note the 
relatively minor differences between the sensitivities, the Feldman-Cousins limits 
and the Rolke-Lopez limits. Our final result is the Rolke-Lopez 90% confidence up- 
per limit including the systematic error shown in the fifth column of the table. The 
single cut check result, which also includes the systematic error, shown in the last 
column of the table, agrees with our final result as well. All modes are (xlO~^). 

Decay Sensitivity F-C R-L R-L incl. Single Cut 





Mode 










ar 


incl. ar 


D+ - 


K+^i~ 




7.5 


11 


9.1 


9.2 


12 


D+ - 


K-fi- 




4.8 


13 


13 


13 


12 


D+ 


vr+yU^ 




7.6 


9.3 


8.7 


8.8 


7.4 


D+ 


vr^/i^ 




5.5 


4.6 


4.8 


4.8 


5.1 


Dt- 


K+iJ.- 




33 


31 


33 


36 


38 


Dt- 


K-fi- 




21 


11 


13 


13 


20 


Dt 


7r+/i^ 




31 


20 


24 


26 


18 


Dt 


—I- TT" 




23 


29 


26 


29 


22 
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Table 4 

FOCUS results compared to other experiments and recent theory. The previous 
limits, except for the E687 K~ fi~^fi~^ [4] are from Fermilab experiment E791 

[5]. The theory estimates come from [2] (SM-1), [3] (SM-3), and [1] (R-Parity MSSM 
and SM-2). Note that the SM estimates from [1] use a formalism close to [3], and 
at present there is some discrepancy in the invariant Mn mass behavior for the SM 
estimates in SM-3 [1] and SM-1 [2]. All modes shown are (xlO"*^). 



Decay 


This 


SM-1 


SM-2 


SM-3 


MSSM 


Previous 


Mode 


Analysis 








R-Parity 
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D+ 


9.2 


0.007 
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D+ K-H+H+ 
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D+ 7r+/x-/x+ 
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1.0 
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1.8 


15 


15 


D+ 7r^/i+/i+ 


4.8 
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D+ ^ K+f,-fi+ 


36 


0.043 
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D+ ^ 


13 










180 


Df ^ 7r+/x-/x+ 


26 


6.1 
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D+ ^ 7r-/i+/i+ 


29 










82 
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